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ABSTRACT 


In this work, the influences of various environmental scenarios on the bottom 
interface and volume reverberation in shallow water were numerically analyzed. Based 
on similar modeling reverberation geometry defined in previous works, the numerical 
analyses were conducted for broadband pulse signals to generate complex reverberation 
structures in the time-domain. The reverberation model used is based on the well- 
documented Parabolic Equation (PE) approximation. The environmental scenarios are 
divided into three main categories. They include different sound speed profiles, different 
levels of bottom interface roughness and different bottom volume fluctuations. While 
one category is being analyzed, the other two are held constant. The various analyses 
include broadband two-way reverberation levels comparisons, vertical correlation 
analysis and power spectral analysis. 
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I. INTRODUCTION 


As the Navy expands its operations towards the littoral region, the ability to 
accurately characterize reverberation in shallow water becomes extremely important. 
Shallow water reverberation affects active sonar systems, impacting the performance of 
underwater detection and tracking systems, as well as acoustic communications. The 
primary mechanisms creating shallow water acoustic reverberation are the propagation, 
the bottom interface and sub-bottom fluctuations, and the rough sea surface. 

In the early 1990’s, the Office of Naval Research (ONR) sponsored a multi-year 
reverberation program known as the Acoustic Reverberation Special Research Project 
(ARSRP). The main goal was to study the primary causes and nature of acoustic 
reverberation in the deep ocean. The area examined was near the Mid-Atlantic Ridge, an 
area where the bottom is covered mostly by highly variable topography and hard rock 
stmctures. The results from the experiments showed that the predominant mechanism 
for high reverberation levels was the interaction of the propagation with the bottom 
interface topography, with harder rocks generating higher scattering levels than softer 
sediment ponds. To a great extent, the general stmcture of the reverberation coincided 
well with the predicted two-way transmission loss, with high reverberation returns 
occurring at areas where the acoustic energy interacted strongly with the bottom. This 
implied that much of the long wavelength statistics of the acoustic energy could be 
predicted by using propagation modeling. 

ONR has recently sponsored another reverberation study within a program 
codenamed ASIAEX, to record reverberation signals and collect oceanographic data in 
the shallow waters of East China Sea. In contrast to the deep ocean reverberation 
experiment, the seabed of this littoral region has a much smoother, softer and more 
penetrable bottom layer such as sand or mud, covering a harder sub-bottom layer of 
coarser sand, gravels or rocks. 

For the past few years, thesis work by previous students^^^’^^^ was focused on 
examining the shallow water bottom reverberation using the Monterey-Miami Parabolic 
Equation (MMPE) propagation model.The MMPE model was first developed by 
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Smith and Tappert in 1994. Since then, the MMPE model has been further improved to 
include bottom interface and volume perturbations in order to create more reahstic 
environmental models. The interface perturbation is simply a spatial displacement 
perturbation while the volume perturbation involves both sound speed and density 
fluctuations within the sediment. 

The previous work modeled a 16-element vertical line array (VLA) geometry to 
support monostatic and bistatic (vertical separation) reverberation computation in a 
shallow water environment. The element located at 48m was chosen as the source with 
all 16 elements acting as receiver. The previous environment was a 100m isospeed water 
column with a single water/bottom interface. The bottom was modeled with an interface 
root-mean-square (rms) roughness of Im and a 15m/s rms volume sound speed 
perturbation. With the geometry and environmental parameters defined, several 
reverberation analyses were then conducted to compare perturbed data with and without 
density fluctuations in the environment for both continuous wave (CW) and broadband 
signals. The analyses included reverberation pressure level comparisons, vertical 
correlation analysis, peak correlation analysis and spectral analysis. Some of the 
significant results were as follows: 

• The influence of volume density fluctuations was to reduce later reverberation 
levels relative to earher levels but did not affect the stmcture significantly. 
This was due to the direct correlation between volume sound speed and 
density fluctuations used in the model. 

• It was noted that the CW analysis was unable to capture coherent stmcture of 

the volume reverberation pressure level due to the inabdity of CW to resolve 
multi-path influence. Therefore, the peak vertical correlation analysis was 

vahd only for broadband pulse computations. 

• The peak vertical correlation analysis suggested that the volume reverberation 

decorrelated across the vertical array more rapidly than interface 

reverberation. This was presumably due to multi-point/multi-depth scatter 
contributions of the volume producing more vertical stmcture than the 
interface. 
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• Spectral analysis of both CW and broadband pulse calculations suggested that 
response of the interface reverberation has a slope on the order of -0.125 for 

both CW and broadband data. However, the volume response showed a 

steeper -0.75 slope for CW and -0.25 slope for broadband signals. 

Taking into consideration the findings obtained from previous work, the objective 
of this thesis is to examine the influence of various environmental profiles on the 
character of the predicted bottom reverberation for broadband signals. The MMPE 
model was used to generate the bottom reverberation data, using the same 16-element 
VLA reverberation geometry defined previously. The MMPE model application program 
used is named “MMPEREVERBDENS2” and is written in the FOTRAN programming 
language. The environmental models used include different sound speed profiles and two 
bottom interfaces (water/bottom interface and bottom/sub-bottom interface) with 
different interface roughness and volume perturbation. The sound speed profiles were 
based on data taken during the recent East China Sea portion of ASIAEX, while the 

bottom and sub-bottom characterization were based on preliminary reports of the geo¬ 

acoustics of the region^"^^. Signal processing and analysis in the time domain were 
performed using MATLAB. Data comparisons were made then with reference to 
reverberation results generated using an environmental profile typical of a shallow water 
region. 
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n. NUMERICAL METHODS AND IMPLEMENTATION 


A. REVERBERATION THEORY 

Reverberation arises from the scattering or re-radiation of the transmitted signals 
from unwanted targets such as marine hfe, bubbles, the sea surface, sea bottom and sub- 
bottom fluctuations. The focus in this thesis is to examine the reverberation generated 
by the water/bottom interface, bottom/sub-bottom interface and the bottom volume 
fluctuations in shallow water. The sea bottom interface variabihty ranges from small 
features producing Bragg scatter to large features such as sea mounts and pinnacles 
producing mostly specular reflection of larger wavelengths. 

The theoretical principle of the scattering mechanism is the same for both 
monostatic or bistatic reverberation. The theoretical treatment will only focus on the 
monostatic mode since its numerical implementation is simpler. 

In decibel units relative to 1 pPa and reference length scale of 1 m, we define the 
mean reverberation pressure level, RPL, as 


RPL = lOlog 


/ 

\l 




- 

— SL + DIj + DIj^ + lOlog 

b 

\Pr^ 

/j 

_ ^0 _ 


-RL, 


'b,v ■ 


( 2 . 1 ) 


where SL is the source level, DIj is the directivity index for the transmitter, D/^ is the 
directivity index for the receiver, A/l^is the ensonified area in the horizontal direction, 
/?Qis the reference distance, and RL,^^ is the reverberation loss per unit area for either the 
bottom interface or the volume. 
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1. Bottom Interface Scattering 

The reverberation loss for the bottom interface is defined by^^^ 


RL,^2TL,-S, 


( 2 . 2 ) 


and 


5,=101og 




3271 - 


wA^K) 


y-V J 


(2.3) 


where TL^ is the average transmission loss (based on long wavelength components) from 
source to the scattering patch at the bottom, is the fuU-wave scattering strength due to 
the small-scale interface roughness, is the wavenumber, is a reference sound speed, 
Ac is the difference in sound speed between the water and the bottom at the interface, and 
W^Q-k^) is the two-dimensional (2-D) spectrum of the interface roughness evaluated at 
the Bragg wavenumber for monostatic reverberation. 


2. Volume Scattering 

The volume reverberation loss cannot be expressed simply in terms of the two- 
way transmission loss but instead must be the integral over depth of the quantity 

at each range r, where n(r,z) is the approximate refractive index 

based on only long wavelength perturbation, and \|/(r,z) is the field function of the two- 

way propagation (as defined in Eq 2.11). The reverberation loss for the volume is then 
defined by^^^ 


and 



(2.4) 


5,=101og 


3271 




(2.5) 
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where is the volume scattering strength, W2,{2kQ) is the 2-D horizontal spectrum of 
the volume fluctuations, which is assumed horizontally isotropic and independent of 
depth evaluated at the Bragg wavenumber. 


B. MONTEREY-MIAMI PARAB OLIC EQUATION (MMPE) MODEL 

The parabohc equation (PE) method is a popular numerical approach for solving 
the acoustic wave equation. The MMPE ModeP^ is based upon the parabolic 
approximation of the wave equation and, therefore, a brief description of this approach 
would be useful. 

We start by representing the time harmonic acoustic pressure field defined in a 
cylindrical coordinate system and assuming azimuthal S 3 mmietry, 

P(r,z,cot)= p(r,z)e”'“'. (2.6) 


Cyhndrical coordinate is chosen because the shallow water sea can be portrayed as a thin 
waveguide on the surface of the earth. Azimuthal symmetry is assumed because the 
ocean environment tends to exhibit weak azimuthal dependence. Substituting Eq. (2.6) 
into the wave equation in cylindrical coordinates leads to the Helmholtz equation, 

V^p(r,z)+ ^ pir,z) = 0, (2.7) 

c{r,z) 

where 


r—+— 
r dr dr dz^ 


( 2 . 8 ) 


The Helmholtz equation can be factored by introducing the operator notation 

ap=(p+e + l)^, (2.9) 

where 

r 1 d^ 

e = n^-l, n = ^ and|i=——, (2.10) 

c dz 


and Co is the reference sound speed typical of the ocean volume. 
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Taking into consideration the effect of cylindrical spreading, and proper factorization of 
the Helmholtz equation, the acoustic pressure may then be defined as. 


p{r,z)^P„ 



( 2 . 11 ) 


where is the envelope function or PE field function. The parabolic equation for the 
field function is then defined by 


dr 




0^0 


( 2 . 12 ) 


where 


( 2 . 13 ) 

is a Ha milt onian- lik e operator which defines the evolution of the PE field function in 
range. 

The relationship between values of \|/ at different ranges can be defined by 

\|/(r-l-Ar) = <l>(r)\|/(r), (2.14) 


where ^>(r) is a propagator that marches the solution out in range. The MMPE model 
employs a split-step Fourier (PE/SSF) method'^^^ to provide a representation of the 
propagator ^>(r). This method is utilized primarily because of the speed and simplicity of 
the PE/SSF method. With the Fast Fourier Transform (FFT), the PE/SSF implementation 
can be represented by 


\\f{r + Ar,z)-e 
where 





e \|/(r,z) 

[ 

V JJ 


(2.15) 


and 


(2.16) 


T =1- 

op 


1 - 


Z 


K o ) 


ni/2 


(2.17) 
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C IMPLEMENTATION OE THE REVERBERATION PROBLEM IN MMPE 

Having described the concepts of reverberation theory and the MMPE model, we 
will now focus on the theoretical treatments for generating perturbation to both interface 
roughness and volume sound speed, including the influence of density fluctuations of the 
volume. The incorporation of these effects into the MMPE model is also discussed. The 
theoretical basis for modeling the interface roughness is based on the work of Goff and 
Jordon,and the development for the volume perturbation theory is based on 
Yamamoto’s work.'^’^ 

1. Interface Roughness 

We assume a two-dimensional (2-D) interface spectrum of the form^^^ 

^ and k, = + iJ , (2.18) 

where kr is the horizontal spatial wavenumber vector, p is a normalization factor, Lcorr is a 
correlation length scale, |3 is the spectral exponent and K and L are the horizontal 
wavenumbers in the x- and y-directions, respectively. 

If the 2-D spectrum W 2 (kr) is assumed to be independent of direction (isotropic), 
then the norm aliz ation factor p can be defined in terms of the root-mean-square (rms) 
roughness a ^ by requiring 

2n]w^{k^)k4k^^G\ (2.19) 

0 

which leads to 

We simply evaluate W 2 (kr) at kr=2ko for the scattering amphtude since it is caused 
by Bragg scatter (evaluated along the line of propagation for monostatic reverberation). 
However, for the long-wavelength interface roughness, we need the full spectrum, i.e. the 
one-dimensional (1-D) spectrum along the x-axis. This can be done by taking the 1-D 
transform of W 2 (K,L) along a slice at y=0, defined by 
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W,{K)^ \w^{K,L)dL. 


( 2 . 21 ) 


In cylindrical coordinates, we have 

k 


i / , (K)dK=r^ (1+iL'f') 


Mi 

2 2 


( 2 . 22 ) 


where 


Y = 



r 

rn 

r 


jz j 


UJ 


2 2 
Ljz_ ±J. 


Ttr 




(2.23) 
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In order to generate a 1-D roughness realization from Eq. (2.21) or Eq. (2.22), we 
transform the 1-D amphtude spectrum that has been scaled by a random amphtude and 
phase. That means we can define the roughness reali z ation as 


ri(x) = 1“ Si(K)e‘^dK, 


where 


S,(K) = [W,(K)]yM{K)e‘^ 


,B(K} 


(2.24) 


(2.25) 


and A and 6 are random numbers for aU values of K. The random phase and amphtude of 
each component can be obtained from Eq. (2.26) and Eq. (2.27), respectively 

e^lKr, (2.26) 


A=^-\n{r2) 


(2.27) 


where both ri and r 2 are independent uniformly distributed random variables in the 
interval [0,1]. In practice, we simply use 


and rescale the result by its rms value, i.e. 


W,(K) = {l + I:„,X) 

r\ix) 


P+i 
2 2 


(2.28) 


(fi"(^)) 
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2. Volume Sound Speed Fluctuations 

The sediment volume sound speed perturbation may be modelled by a three- 
dimensional (3-D) volume spectrum given by*^'^ 

W,{K,L,M) = A" (k^+ + (2.29) 

27t 

where B is the spectral strength constant, p is the spectral exponent, and A = — = — is 

the horizontal-to-vertical aspect ratio describing the anisotropy of fluctuations in the 
sediment, K and L are the horizontal wavenumbers in the x- and y-directions, 
respectively, and M is the vertical wavenumber. 


To evaluate the reverberation due to the volume perturbation, we need an 
expression for the 2-D horizontal spectrum (assuming strongest scattering near 
horizontal). It is defined as 

IT 2 (/^, L) = IIT 3 {K, L, M)dM . (2.30) 

Substituting Eq.(2.29) into Eq.(2.30), we have 

= + dM . (2.31) 

71 Q 

For P = 2, Eq. (2.31) can be reduced to 

+ . (2.32) 

For the values of 5-5x10“^ and A~5 chosen from Yamamoto’s findings,^’^ Eq. (2.32) 
reduces to 

_ 3 

W^{K,L)=a[K^ + L^\ (2.33) 

where 
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(2.34) 


a = — = 5x10 ^ 

2A 

However, for the forward propagation, we need only the 2-D vertical spectrum in 
the {r,z) plane. It can be defined by 

W'^{K,M)^ \%{K,L,M)dL^^ - f Ta" (/f" + Z?) + M"T " dL. (2.35) 

For p = 2 , Eq. (2.35) becomes 


W^(K,M) =a'[25Z:" + M"] “ 2 , 


(2.36) 


where 


a'= 1.25x10“' . (2.37) 

To generate 2-D vertical volume sound speed fluctuation realizations, we define a 
reaUzation as 

Cpb (x,z) = jj S^(K ,, (2.38) 

where 

S^{K,M) = [W'(K,M)f^A(K,M)d^^^’^’ . (2.39) 

Notice that since we have treated the sound speed perturbation in the volume in the 
vertical, we are really generating a series of vertical realizations at each range step. 

In similar fashion to the interface, the 2-D random phase and the amphtude 
variations can be obtained, respectively, by 

%{K,M)^llir,{K,M) (2.40) 

and 

AiK,M) = yl-\nir,(K,M)) (2.41) 

where both r^{K,M) and r 2 {K,M) are now a matrix of uniformly distributed random 
numbers in [0,1]. In practice, we use 


12 



(2.42) 


_ i 

+M^) 2 

and rescale by the appropriate rms values. 

Eqs. (2.28) and (2.42) are the results for the interface roughness and volume 
perturbation, respectively. These are the generic spectral models used in generating the 
reahzations for implementation in the MMPE model. 


3. Density Fluctuations in Sediment 

Variability in density, p, is incorporated into the PE model by defining the 
effective index of refractiorf^^ 


/2 2 , 
n -n +■ 


2kl 


1 2 3 


2 

-V^p-- 

-Vp 


P 2 

P 



(2.43) 


Consistent with the numerical treatment that assumes the environment is range- 
independent over a range step, and the fact that sediment properties are largely 
horizontally stratified, we may simplify this to 


/2 2 , 
n —n +■ 


2e 


1 a^p 

p dz^ 




1 dp 

vP 


(2.44) 


For the forward problem, the sound speed index of refraction is based only on large scale 
features, such that 


and 


2 ^ 2 \ 

{r,z) 



(2.45) 


(! + (?, +5,) = c +5c, , 5c, =c, (h^+5,) , 


(2.46) 


where c,^ is the mean bottom sound speed at the interface, - y is the normalized 

gradient of bottom sound speed, and 5, is the zero-mean random perturbation for the 
long wavelength component. 
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According to the analysis of Yamamoto (1996)/^^ the relative fluctuations in 
density relate to the relative fluctuations in sound speed according to 


Sp _ 2(p,-po) 5c 

Po ^Pq —Pr Cq CJ) 


(2.47) 


where 


CQ=Cu{\+b^) andy = 


P.-Po 

2po-p, 


(2.48) 


Notice that Pgand Cq are the averaged values of the density and sound speed in the 
sediment and = 2650^^/3 is the density of the grain. We may then write 


p = Po+5p = Po 1 + —5c 


(2.49) 


Taking the first and second partial derivatives of Eq. (2.49) with respect to depth, z, and 


neglecting depth gradients in either c^or pQ , we obtain 

dz Cq dz 


(2.50) 




(2.51) 


The sound speed fluctuation is defined by 


where 


5 c^ = CqS (x, z) = jj52 (/s:,M y^^e‘^^dKdM , 
S, {K,M) = [w;{K,M)Y^ A{K,M) . 


(2.52) 


(2.53) 


Substituting Eq. (2.52) into Eq. (2.50) and Eq. (2.51) yields 


^ ff [iM )S2 {K, M)e‘''^e‘^^dKdM 

dz Cn 


(2.54) 
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^ ff i-M^)S 2 ( K, M)e‘^e‘^^dKdM . (2.55) 

dz Cg 


The sediment effective index of refraction can now be derived by substituting Eq. (2.54) 
and Eq. (2.55) into Eq. (2.44), which becomes 

nl^{x,z)^n^ {x,z) + -^\cl{x,z)+ P(x,z)1, (2.56) 

2^0 

where 

a (X, z) = , \\m^S^[K, M) (2.57) 

Co(x,z)-'-’ 


and 


P (x,z) = + ^ 


2y 


Co(x,z) 


|j MS^ {K, M)e‘^e‘^^dKdM 


(2.58) 


These parameters will be computed in parallel with 6 in the MMPE model. 
In the sediment, there is then the additional propagator term 

where 


(2.59) 


Up {x,z) = --^[oc (x, z) + p (x, z)] . 

4/Cn 


(2.60) 
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D. TIME-DOMAIN PROCESSING 

The theoretical treatment of the reverberation loss for the bottom interface RLy 
and the volume RLv discussed previously were based on CW analysis. However, our 
focus here is on broadband analysis. To predict the effect of pulse propagation in time, 
we need to run the MMPE model over a spectrum of frequencies. This allows us to treat 
coherent interference effects by separating the multipath effects. The time-domain 
analysis of the interface and the volume reverberation will provide the general picture of 
the two-way travel time stmcture of the reverberation loss. We can then determine the 
reverberant field at each range step and continue the propagation through the entire water 
column of interest. 

1. Time-Domain Analysis of the Interface 



Figure 1. Two-Way Return from a Scattering Patch 


The geometry of a two-way return from a scattering patch adapted from Smith 
and Cushman'^^^ is shown in Figure 1. The travel time structure of the one-way forward 

r 

propagating field at range step is denoted by p+(rm,z,t) where t = T - — is the reduced 

Co 

time and T is the actual travel time. The two-way pressure field at the receiver is the 
convolution of two, one-way fields in the tim e-domain^^^ 

= (2.61) 

where and are the forward propagated pressure fields from the transmitter and 
receiver to the scattering point evaluated at the bottom interface, respectively. Note that 
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here x = T - — is the reduced time of the reverberation and T is the actual travel time. 

Furthermore, the receiver and transmitter need not be co-located in the water column. By 
reciprocity, the propagated field from the receiver to the scattering point, is the 

same as the propagated field from the scattering point to the receiver . 


The time-domain convolution of the two field functions is also the scalar 
multiplication of these functions in the frequency domain. Hence, the two-way field in 
the frequency domain from the interface can be expressed as 

P2-«ayi,ir„J) = P^n(r„,f)p +Rb /), (2.62) 

where 


P^Tbir^J) 




(2.63) 


and 


P+Rb^^m’f^ I - 




+Rb^ m 


(2.64) 


The two-way travel time structure of the reverberation loss for the bottom interface, RLt, , 
due to a single bottom patch can then be defined as 

P-.fr.,.') = Aj ,/)«-""# , (2.65) 


where the constant A is included to incorporate aU the other factors needed to define 
reverberation loss, RLb. This provides the two-way travel time stmcture due to scattering 
from range r^. The calculation is then continued for each range step. The total field at the 
receiver is computed by coherently summing up all the pressure values from the different 
range segments, r^, by matching up the discrete arrival times, according to 

M 

P-b{K)='^P-b{^mAn)^ ( 2 - 66 ) 

m=l 

where p.b is the total interface reverberation pressure received at the receiver at time 
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2. Time-Domain Analysis of the Volume 

The reverberant field due to each depth/range point is computed by combining the 
source-to-patch and patch-to-receiver propagating field, according to 

P2-way.v(^m^z ,/) = n(r^, z)p^j{r^,z J)p^R{r^Z ,/), (2.67) 

where the two-way reverberation signal is computed for every grid point of interest 
(z> Zj, always) at a particular frequency, /, and n{r^,z) is the local index of refraction at 

the grid point. The reason for multiplying by the local index of reffacttion, n(r^,z), is to 
provide the same weighting used in the CW treatment. 

Fourier transform gives the tim e-domain response 

p, (r,z,t) = fifp (r ,z,f)e~‘^^^df. (2.68) 

The two-way travel time sttucture of the reverberation loss for the volume, RLy, can then 
be derived from 

j P2-wayA^m^Z,t)dz , (2.69) 

Z>Zb 

which is the coherent sum over all depths below the interface at range step m. Note that 
the constant B is included to account for all the other terms needed to define 
reverberation loss, RLy. The single set of time series can then be matched and summed to 
give 

M 

= ( 2 - 70 ) 

m=l 

where p.y is the two-way time domain pressure defining the volume reverberation loss at 
the receiver at time, due to the entire volume of interest. 
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m. MODELING GEOMETRY & ENVIRONMENT 


A. MULTI-STATIC REVERBERATION GEOMETRY 

Since the varying spatial properties of the broadband signal are of interest, a 
vertical line array (VLA) with 16 elements was chosen. The VLA will support both 
monostatic and bistatic (in the vertical) reverberation computations/measurements. The 
array was located vertically in the lOOm deep water column and each element of the VLA 
was assumed to be a point source/receiver. Figure 2 shows the VLA and its geometry 
relative to a scattering patch at a horizontal range, r. From Figure 2, we see that the 16- 
element VLA spans the water column from 20m to 80m with 4m separation in depth 
between adjacent elements. A single element located at 48m depth was chosen as the 
source with all 16 elements receiving the reverberation. To perform broadband analysis. 



Figure 2. Geometry of VLA and Scattering Patch 
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B. THE ENVIRONMENTAL MODELS 


The maximum propagation range is 5km. The mean bottom depth is lOOm while 
the mean sub-bottom or deep-bottom depth is 105m. Hence, there are two bottom 
interfaces, the water/bottom and bottom/sub-bottom (or bottom/deep-bottom) interfaces. 
This is physically representative of a smoother, softer and more penetrable bottom layer 
such as sand or mud, covering a harder sub-bottom layer of coarser sand, gravel or rocks. 

Seven sets of environmental parameters are used to run the MMPE model. The 
main differences between each set are the variations in sound speed profiles (SSP) of the 
water column, water/bottom and bottom/sub-bottom interface roughness and bottom/sub¬ 
bottom volume sound speed perturbation. One of the seven profiles is a reference model 
that follows the typical profiles of what is expected in the shallow water region. The six 
other models are divided into three categories based on the differences mentioned. 
Comparison of results within each category will be done with respect to the reference 
model. 

1. Variations in Sound Speed Profiles 

Six representative sound speed profiles were obtained from the ASIAEX 
experiments and are plotted in the figure below. Of the six, the two extreme left and right 
sound speed profiles (SSPl and SSP2) and the average sound speed profile (SSPAvg) 
were chosen for the modeling runs. The SSPAvg is used in the reference model. 



Figure 3. Typical Sound Speed Profiles of East China Sea 
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2. Variations in Interface Roughness 

Three sets of interface root-mean-square (rms) roughness are used. They are 
(0.5m, Im), (Im, 2m) and (2m, 4m) for each water/bottom and bottom/sub-bottom 
interface pair. Pbts of the rough interface realizations for the bottom and sub-bottom 
bathymetry for different rms values are shown in Figure 4. 



RsnO@ r^:Tri? 


Inlorfflce 11115 RDughnpaa 




(kin) 


Figure 4. Plots of Different Interface Roughness 
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From Figure 4, increased roughness is observed as the rms values are increased. 
The longer wavelength perturbation also become more prominent with increased r ms 
roughness. It is also observed that the largest rms roughness values (2m, 4m pair) 
produce some sub-bottom protmsions through the water/bottom interface. It should be 
noted that the vertical depth scale is in meters and the horizontal range scale is in 
kilometers, which makes the plots look extremely rough. The (Im, 2m) pair is used in 
the reference model. 

3. Variations in Volume Perturbations 

In order to depict a softer bottom layer covering a harder sub-bottom layer, the 
sound speed used in the bottom layer is IVOOm/s while that of the sub-bottom layer is 
1760m/s. This is typical of a sand/mud layer covering a gravel/rock sub-layer. The rms 
volume sound speed perturbations chosen are 5m/s, 15m/s and 45m/s. The 15m/s rms 
sound speed perturbation is the reference parameter. 

The effects of the volumetric rms perturbation to the sound speed are illustrated in 
the figures on the next page. Apparent from the plots is the increasing contrast in the 
sound speed profile for an increasing volume sound speed perturbation. The rms sound 
speed perturbation of 15m/s is used in the reference model. It is not likely that the 
volume sound speed profile would vary considerably due to inhomogeneities and mixture 
of different sediment types as in the case with rms sound speed perturbation of 45m/s. 
However, the rms sound speed perturbation of 45m/s is still used for comparison 
purposes. 
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Figure 5. Plots of Sound Speed Data with Variations in Volume Perturbations 
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The differences in each environmental profile are summarized in the Table 1. 


Environmental 

Profiles 

Sound 

Speed 

Profiles 

Interface Roughness (rms) 

Bottom Volume 
Sound Speed 
Perturbation 
(rms) 

Category 

Water/ 

Bottom 

Bottom/ 

Sub-Bottom 

1 

SSPAvg 

Im 

2m 

15m/s 

Reference 

Model 

2 

SSPl 

Im 

2m 

15m/s 

Variations in Sound 
Speed 

3 

SSP2 

Im 

2m 

15m/s 

Variations in Sound 
Speed 

4 

SSPAvg 

0.5m 

Im 

15m/s 

Variations in Interface 
Roughness 

5 

SSPAvg 

2m 

4m 

15m/s 

Variations in Interface 
Roughness 

6 

SSPAvg 

Im 

2m 

5m/s 

Variations in Volume 

Perturbations 

7 

SSPAvg 

Im 

2m 

45m/s 

Variations in Volume 
Perturbations 

Table 1. Environmental] 

Profiles 


4. Other Parameters 

The other parameters used in execution of the MMPE model, with their respective 


resident input filenames are specified below: 


Filename/Parameter 

Value 

Remarks 

Main Control File: pefiles.inp 



Number of depth points 

256 

Radix-2 integer required for FFT 

Minimum depth 

0 m 


Maximum depth 

400 m 


Number of range steps 

833 


Minimum range 

0 m 


Maximum range 

5.0 km 


Range step size 

6 m 


Maximum computed depth 

400 m 


Reference sound speed 

1500 m/s 


Source File: pesrc.inp 



Source depths 

Varying 

Array elements at 20, 24, 28, 32, 36, 

40, 44, 48, 52, 56, 60, 64, 68, 72, 

76 and 80 m depths. 

Center frequency 

250 Hz 


Frequency bandwidth 

250 Hz 


No. of Frequencies 

512 

Radix-2 integer required for FFT 

Sound Speed File: pessp.inp 



Water column sound speed 

3 sets 

SSPl, SSP2 and SSPAvg 

Range independent 

No. of SSPs points 

58 


Bathymetry: pebath.inp 
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Filename/Parameter 

Value 

Remarks 

Mean bottom depth 

100 m 

Range independent 

No. of depth points 

1 


Bottom properties: pebotprop.inp 



Bottom sound speed 

1700 m/s 


Sound speed gradient 

1 /s 


Relative density 

1.6 

No density variation 

Compressional attenuation 

0.15 dB/km/Hz 


Shear speed 

0 

Negligible 

Shear attenuation 

0 

Negligible 

Sub-Bottom Bathymetry: pedbath.inp 



Depth 

105 m 


Sub-Bottom Properties: pedbotprop.inp 



Sub-bottom sound speed 

1760 m/s 


Sound speed gradient 

1 


Relative density 

2 

No density variation 

Compressional attenuation 

0.2 dB/km/Hz 


Shear speed 

150 


Shear attenuation 

0.5 


RMS Perturbations (input during running of model) 



Water/bottom interface roughness 

3 sets 

0.5m, Im, 2m 

Bottom/sub-bottom interface roughness 

3 sets 

Im, 2m, 4m 

Volume sound speed fluctuation 

3 sets 

5m/s, 15m/s 45m/s 


Table 2. Input Environmental Parameters 
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IV. POST PROCESSING AND RESULTS 


A. POST PROCESSING 1 - TRANSMISSION LOSS 

Before examining the two-way reverberation results, it is useful to examine the 
one-way transmission loss plots for the various environmental profiles to provide a 
preliminary examination of what may be expected. The transmission loss plots were 
based on propagation of CW signal with frequency of 250Hz from a source at 48m depth. 


1. Transmission Loss Due to Variations in Sound Speed Profiles 

The transmission loss plots for the three different sound speed profiles were very 
similar in stmcture. Hence, there should be httle variation in the reverberation loss. It is 
noted that energy is being refracted upwards in the sub-bottom layer. This is due to the 
sound speed gradient of Im/s/m included in the sub-bottom layer. 



Figure 6. Transmission Loss Due to Variations in Sound Speed Profiles 


2. Transmission Loss Due to Variations in Interface Roughness 

The transmission loss stmctures for the three different sets of interface roughness 
look similar, as seen in Figure 7. However, careful observation shows that the higher 
interface rms roughness values have lower transmission loss. This may be because the 
higher roughness values have caused the sediment layer between the water/bottom and 
the bottom/sub-bottom interfaces to be thinner and allow the acoustic energy to interact 
more readily with the denser and faster volume below the bottom/sub-bottom interface. 
Thus, more forward propagation of the energy occurs. In the case of the largest 
roughness values, the sub-bottom protmsions through the water/bottom interface may 
have cause even greater forward propagation of the energy. Therefore, lower 
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reverberation loss (or higher reverberation level) for higher interface rms roughness 
values may be expected. 



Figure 7. Transmission Loss Due to Variations in Interface Roughness 

3. Transmission Loss Due to Variations in Volume Perturbations 


It is observed that the transmission loss is greatly affected by the variations in 
volume perturbations. The high sound speed fluctuation value of 45m/s increases the 
transmission loss significantly. The transmission loss from the 15m/s volume sound 
speed fluctuation is also considerably higher than that of the 5m/s fluctuation. In an 
attempt to determine the cause, the sound speed data plot from Figure 5 is examined. It is 
found that at short ranges there is a patch of slower bottom sound speed between 0 and 
0.3km. This patch of slower bottom sound speed has allowed more of the signal to 
penetrate the bottom. 



Figure 8. Transmission Loss Due to Variations in Volume Perturbations 


In order to further examine the high transmission loss created by the high volume 
sound speed fluctuation, the model is run with a new random seed to create a different 
bottom sound speed profile (Figure 9) with a sound speed fluctuation of 45m/s. It is 
found that at short ranges (0 to 0.3km), the bottom sound speed profile is shghtly higher 
than the previous. The transmission loss plot also shows lower transmission loss out to 
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that range. Note that there is another slow patch further out (about 0.6km to 1km) that 
causes high transmission loss at a further range. This confirms that a high volume sound 
speed fluctuation may create a slow bottom that significantly reduces reflection, thus 
increasing transmission loss. We should then expect to see high reverberation loss at 
high volume sound speed fluctuation. It should be noted that the rest of the results and 
analysis in the following sections are based on the original sound speed profile data 
created by the original seed. 


Figure 9. 
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Sound Speed Profile Data and Transmission Loss for New Random Seed 


B. POST PROCESSING 2 - REVERBERATION LOSS 


1. Time Domain Reverberation Analysis 

To evaluate the bottom reverberation loss, analyses of the pulse propagation in the 
time-domain for both the reverberation stmcture of the interface and volume were 
necessary. The time-domain analyses were previously explained in Chapter E, Section D. 
The equations required to formulate the MATLAB implementation for the interface and 
the volume reverberation loss are summarized below: 

• For the interface reverberation loss, RL^, 



Pl-.ayJ, ’ /) = P^Tb (fm > f)P^Rb /) » 

(4.1) 



(4.2) 


M 

P-b{K) = T.P-b{rm,tn)^ 

m=l 

(4.3) 


«Z.,((.)=-201og[|p_,((.)|]. 

(4.4) 
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• For the volume reverberation loss, RL^, 



,/) = n(r„, z)p^Tir^,z,f)p^„ir^,z ,/), 

(4.5) 


p, (r , z,t)- [ p, (r ,z ,f )e~'^^'df , 

^ 2-way, V ^ m’ ’ ^ J ^ 2-way,v ^ m ’ ^ ^ ^ 

(4.6) 


P_,(r„,0= j P2-^ayA^m^Z^^)dZ , 

(4.7) 


Z>Zt 



M 



p-v(l) = Ef-v((;>4)’ 

m=l 

(4.8) 


M(r.)=-201og[|p_,(*.)|]. 

(4.9) 


2. Reference Model Reverberation Loss 

With a source depth of 48m, the general color maps of the water/bottom and 
bottom/sub-bottom interfaces and volume reverberation losses (for 2-way transmission) 
were obtained from the reference model pgure 10). Fine multipath structures due to the 
interface and volume fluctuations were present. Careful observation shows that the 
arrival times of the first signal returns occur earber at the receiver at 80m, while the latest 
arrive at the receiver at 20m. This was due to the shorter distance between the deeper 
receiver and the bottom, thus receiving returns earber than the shaUower receivers. The 
color map of the bottom/sub-bottom interface reverberation loss shows a few maxima and 
minima 



Figure 10. Color Maps of Reference Model Reverberation Loss 
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With the receiver chosen at 40m depth, the water/bottom interface reverberation 
and volume reverberation plots are shown in Figure 11. The plots have rather s imil ar 
stmcture. The bottom/sub-bottom interface reverberation has a few conspicuous maxima 
and minima not observed in the other curves. This is due to the sound speed gradient 
included in the sub-bottom layer, which causes the energy to be refracted upwards and 
interact with the bottom/sub-bottom interface. 



Figure 11. Reference Model Reverberation Loss 


3. Reverberation Loss Due to Variations in Sound Speed Profiles 

The reverberation loss color maps for the different sound speed profiles show 
very similar stmctures. This is probably due to the relatively small variations between 
the sound speed profiles. 



Figure 12. Color Maps of Reverberation Loss Due to SSPl 
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Figure 13. Color Maps of Reverberation Loss Due to SSPAvg 



Figure 14. Color Maps of Reverberation Loss Due to SSP2 


With the receiver at 40m depth, the reverberation loss stmctures for the 3 sound 
speed profiles are seen to be quite similar. This is consistent with the color maps. 



Figure 15. Reverberation Loss Due to Variations in Sound Speed Profiles 
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4. Reverberation Loss Due to Variations in Interface Roughness 

The reverberation color maps look similar for the different values of interface rms 
roughness. But when observed carefully, the color maps from the bottom/sub-bottom 
interface reverberation display a more obvious difference, with the higher roughness 
values producing higher reverberation levels. As explained earher, this may be because 
the higher roughness values cause the sediment layer between the water/bottom and the 
bottom/sub-bottom interfaces to be thinner and allow the acoustic energy to interact more 
readily with the denser and faster volume below the bottom/sub-bottom interface. 



Figure 16. Color Maps of Reverberation Loss due to Water/Bottom and Bottom/Sub-Bottom 

Interface rms Roughness of 0.5m and Im Respectively 



Figure 17. Color Maps of Reverberation Loss due to Water/Bottom and Bottom/Sub-Bottom 

Interface rms Roughness of Im and 2m Respectively 



Figure 18. Color Maps of Reverberation Loss due to Water/Bottom and Bottom/Sub-Bottom 

Interface rms Roughness of 2m and 4m Respectively 
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The effect of the variations in interface rms roughness is not very prominent in the 
water/bottom interface and volume reverberation loss, though both showed that the 
higher rms values produce slightly higher reverberation level. However, this effect can 
be seen most clearly from the bottom/sub-bottom interface where the higher interface r ms 
roughness produces reverberation levels which are about lOdB higher. This is probably 
caused by the bottom/sub-bottom interface being nearer the water volume at higher 
roughness values. This allows more energy to interact with the denser and faster sub- 
bottom volume, hence, creating more reverberation. 



Figure 19. Reverberation Loss Due to Variations in Interface Roughness 


5. Reverberation Loss Due to Variations in Volume Perturbations 

It is very obvious from the color maps that an increase in volume rms sound speed 
perturbation reduces the reverberation level for the water/bottom, bottom/sub-bottom 
interfaces and the volume. The reduced reverberation level arising from the increase in 
volume sound speed fluctuations is most probably caused by the patch of slower bottom 
sound speed at short ranges, which allows more energy to penetrate the bottom. This is 
consistent with the transmission loss results. 



Figure 20. Color Maps of Reverberation Loss Due to Volume Perturbation with rms Sound 

Speed Fluctuation of 5m/s 


34 



















































Figure 21. Color Maps of Reverberation Loss Due to Volume Perturbation with rms Sound 

Speed Huctuation of 15m/s 



Figure 22. Color Maps of Reverberation Loss Due to Volume Perturbation with rms Sound 

Speed Huctuation of 45m/s 


The reverberation plots with receiver at 40m depth show that the reverberation 
level is significantly reduced as the volume rms sound speed perturbation is increased, 
which is consistent with the color maps. As explained earher, this is probably the result 
of more energy penetrating the patch of slower bottom sound speed at short ranges. 



Figure 23. Reverberation Loss Due to Variations in Volume Perturbations 
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C POST PROCESSING 3 - VERTICAL CORRELATIONS 

As we are dealing with broadband signals, the vertical correlations for interface 
and volume reverberation were computed in the time-domain. 

1. Vertical Correlation Due to Variations in Sound Speed Profiles 

The color maps of the vertical correlation for the different sound speed profiles 
show very similar structure. It is observed that the reverberation from the bottom/sub¬ 
bottom interface decorrelates fastest, followed by reverberation from the water/bottom 
interface and then the volume. From the color maps, it can also be seen that the deeper 
receivers are receiving signals earlier than the shallower receivers, thus resulting in the 
positive time lag with reference to the source at 48m. 



Figure 24. Vertical Correlation Due to SSPl 



Figure 25. Vertical Correlation Due to SSPAvg 



Figure 26. Vertical Correlation Due to SSP2 
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2. Vertical Correlation Due to Variations in Interface Roughness 

From the color maps, it is noted that the reverberation from the bottom/sub- 
bottom interface decorrelates the fastest across time, followed by the water/bottom 
interface, then the volume, as was observed for the variations in sound speed profiles. 
The various stmctures look rather similar when comparing the various rms roughness 
values with the exception of the bottom/sub-bottom interface, which shows a prominent 
difference. The bottom/sub-bottom interface reverberation displays that an increase in 
the vertical correlation with an increase in the rms roughness. This may be because 
more energy is interacting with the bottom/sub-bottom interface as the rms roughness is 
increased, as previously described. 



Figure 27. Vertical Correlation due to Water/Bottom and Bottom/Sub-Bottom Interface rms 

Roughness of 0.5m and Im Respectively 



Figure 28. Vertical Correlation due to Water/Bottom and Bottom/Sub-Bottom Interface rms 

Roughness of Im and 2m Respectively 



Figure 29. Vertical Correlation due to Water/Bottom and Bottom/Sub-Bottom Interface rms 

Roughness of 2m and 4m Respectively 
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3. Vertical Correlation Due to Variations in Volume Perturbations 

From the color maps, it is observed that the reverberation decorrelates fastest over 
time when the volume sound speed perturbation is increased. This is true for aU 
reverberation types, the water/bottom and bottom/sub-bottom interfaces and the volume. 
It is seen that the reverberation from the bottom/sub-bottom interface decorrelates the 
fastest, followed by the water/bottom interface and then the volume, as before. 



Figure 30. Vertical Correlation Due to Volume Perturbation with rms Sound Speed 

Fluctuation of 5m/s 



Figure 31. Vertical Correlation Due to Volume Perturbation with rms Sound Speed 

Fluctuation of 15m/s 



Figure 32. Vertical Correlation Due to Volume Perturbation with rms Sound Speed 

Fluctuation of 45m/s 
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D. POST PROCESSING 4 - PEAK VERTICAL CORRELATION 


In this section, the peak correlation values were extracted from the vertical 
correlation of the previous section in order to see how the signals decorrelate across 
depth. 


1. Peak Vertical Correlation of Reference Model 

From the figure below, it is seen that volume reverberation decorrelates the fastest 
over depth. The bottom/sub-bottom interface reverberation shows better correlation 
initially but is consistent with the water/bottom interface reverberation when the 
source/receiver separation is increased. 



Figure 33. Peak Vertical Correlation of Reference Model 
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2. Peak Vertical Correlation Due to Variations in Sound Speed Profiles 

The peak vertical correlation structures for various sound speed profiles are very 
similar. Thus, the influence of the water coulum sound speed variability on reverberation 
coherence appears to be minimal. 



Figure 34. Peak Vertical Correlation Due to Variations in Sound Speed Profiles 


3. Peak Vertical Correlation due to Variations in Interface Roughness 

The peak vertical correlation curves for the water/bottom interface due to the 
different interface rms roughness have only slight differences in their stmctures. 
However, for the bottom/sub-bottom interface, the reverberation decorrelates more 
rapidly for increased interface roughness. It seems that the increase in interface 
roughness has caused the vertical coherence of the bottom/sub-bottom interface to be lost 
in the volume. As for the volume reverberation, the decorrelation rate was somewhat 
affected by changes in the interface roughness, but not to the same degree and not 
consistendy. 



Figure 35. Peak Vertical Correlation Due to Variations in Interface Roughness 
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4. Peak Vertical Correlation Due to Variations in Volume Perturbations 

For the case of the water/bottom and bottom/sub-bottom interface reverberation, 
the smaller the volume sound speed perturbation, the faster the decorrelation rate. 
However, the increase in vertical correlation for large volume perturbations is not due to 
a tme increase in structural coherence but rather a significant decrease in reverberation 
levels. In other words, if no signal is received, then it will correlate very weU with itself 
but has no physical significance. For the volume reverberation, the rate of decorrelation 
was not affected as much or as consistently. 



Figure 36. Peak Vertical Correlation Due to Variations in Volume Perturbations 
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E. POST PROCESSING 5 - POWER SPECTRAL DENSITY 

In order to extract the spectral components of the reverberation data for the 
broadband signals, the power spectral density (PSD) and power ratio spectral density 
(PRSD) were computed. The PRSD analysis will be discussed in the next section. The 
magnitude-squared of the ranged-reduced reverberation was analyzed using the Discrete 
Fourier Transform (DFT), such that the PSD is defined by*^*^ 

PSD = DFT[\p_f .\rf], (4.10) 

where p- is the reverberation field of the interface or volume. It should be noted that the 
time-domain pressure field is first converted to range dependent data by the use of the 
reference sound speed co , such that range = time x cq. 

1. PSD of Reference Model 

The PSD plot shows a gradual drop in normalized power in the lower 
wavenumber of up to about O.lm'^ (about 60m length scale) for the interface and volume 
reverberations. After that, a significant drop-off in normalized power was observed. 
The spectral slope over the length scales of 60m to 600m is about -0.55, while the slope 
over the length scales of 10m to 60m is about -1.7. 



Figure 37. Power Spectral Density of Reference Model 
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2. PSD Due to Variations in Sound Speed Profiles 

The PSD structures for the various sound speed profiles are very similar to one 
another for the interface and volume reverberations. Only minor deviation is observed 
for the volume reverberation. Thus, the effects of water column sound speed variations 
on the signal statistics appear to be minimal. 



Figure 38. Power Spectral Density Due to Variations in Sound Speed Profiles 


3. PSD Due to Variations in Interface Roughness 

The water/bottom interface reverberation PSD plot shows little deviation due to 
different levels of interface roughness. The bottom/sub-bottom interface reverberation 
PSD plot shows that the higher interface rms roughness produces slightly higher 
normalized power at wavenumber up to about 0.2m'^ (about 30m length scale), after 
which the lower interface rms roughness contributes to higher normalized power at the 
drop-off. For volume reverberation, the higher interface r ms roughness consistently 
produces lower normalized power throughout. This may be due to the energy interacting 
more with the rougher interface causing more forward propagation, and less energy is 
interacting with the volume. 



Figure 39. Power Spectral Density Due to Variations in Interface Roughness 


43 






































































4. PSD due to Variations in Volume Perturbations 

The water/bottom interface reverberation produces normabzed power which 
decreases with an increase in volume sound speed perturbation, though the 5in/s and 
15m/s rms sound speed perturbation show minor deviation. The stmctures of the 
bottom/sub-bottom interface reverberation PSD plot are similar to that of the 
water/bottom interface. The PSD structures produced by the volume reverberation are 
different. The rms sound speed perturbation of 15m/s produces higher normalized power 
for most wavenumber scales except from 0.02m * to 0.2m * (length scales of about 30m to 
300m). The normalized power from the 45m/s rms sound speed perturbation is 
consistently lower throughout as may be expected from the transmission loss analysis. 



Figure 40. Power Spectral Density due to Variations in Volume Sound Speed Perturbations 
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r. POST PROCESSING 6 - POWER RATIO SPECTRAL DENSITY 

In order to further examine spectral content, the power ratio spectral density 
(PRSD) was implemented. The PRSD is defined in terms of Ihe ratio of the reverberation 


filed of a particular model to that of the unperturbed model, specifically^ 



(4.11) 


where ^perturbed is the reverberation field of the interface or volume of the model of 
interest and ^unperturbed is the reverberation field of the reference model without interface 
roughness and volume fluctuations. 

1. PRSD of Reference Model 

The PRSD of the reference model shows varying results from the water/bottom 
interface, bottom/sub-bottom interface and volume. The bottom/sub-bottom interface 
shows the strongest normalized power, followed by the volume and then the 
water/bottom interface. The normalized power of the bottom/sub-bottom fluctuates 
significandy, followed by the water/bottom interface and then the volume. Over the 
length scales of 10m to 60m, the response of the bottom/sub-bottom interface has a 
steeper slope of about -0.3 compared to a slope of about -0.2 for bottom/sub-bottom 
interface and volume. 


ReTerEfLja MadelPRSD, Sauice. . F?!ec:EU« r. HDm 


— Wacer-EflitBrYi Ifierlace 

- hlEffacE 

- VnlumE 



WavenniiLiEr. jm'^j 


Figure 41. Power Ratio Spectral Density of Reference Model 
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2. PRSD due to Variations in Sound Speed Profiles 

The PRSD structures of the interfaces and volume show little difference due to 
the variations in the sound speed profiles of the water column. As before, this suggests 
that the effect of the different sound speed profiles is minimal. 



Figure 42. Power Ratio Spectral Density Due to Variations in Sound Speed Profiles 


3. PRSD Due to Variations in Interface Roughness 

The PRSD due to different levels of interface roughness show varying results. At 
the water/bottom interface, the PRSDs have about the same normalized power, although 
the fluctuation is greater with the higher interface roughness. At the bottom/sub-bottom 
interface, the higher interface roughness value produces higher normali z ed power, which 
is consistent with the transmi s sion loss result. The volume PRSD shows lesser 
fluctuation with the highest normalized power produced by the interface roughness pairs 
(Im, 2m) followed by (2m, 4m) then lastly (0.5m, Im). 
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Figure 43. Power Ratio Spectral Density due to Variations in Interface Roughness 
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4. PRSD Due to Variations in Volume Perturbations 

The variations in volume perturbations produce significantly different PRSD. At 
the water/bottom interface, the normalized power is highest with the sound speed 
perturbation of 45m/s. The 5m/s and 15m/s sound speed perturbations produce about the 
same magnitude of normalized power, except that there is a significant fluctuation with 
the 5m/s sound speed perturbation. This is not consistent with the transmission loss 
result and the cause is not known at this time. At the bottom/sub-bottom interface, the 
5m/s sound speed perturbation has the highest normal iz ed power, followed by the 15m/s 
and then the 45m/s sound speed perturbations. In this case, the result is consistent with 
the transmission loss analysis, where the loss is greatest at higher volume perturbation. 
In the volume, higher sound speed perturbation produces higher normahzed power, and 
the 5m/s sound speed perturbation shows greater fluctuation than the other two. Again, 
this is not consistent with the transmission loss result. 



Figure 44. Power Ratio Spectral Density Due to Variations in Volume Perturbations 
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V. SUMMARY 


The main focus of this thesis was to investigate the influence of different 
environmental profiles on the reverberation structures due to the water/bottom interface, 
bottom/sub-bottom interface and bottom/sub-bottom volume. Simulations were 
conducted using the same reverberation geometry and environmental parameters as 
defined in previous work. A vertical fine array (VLA) with 16 elements was chosen to 
provide the reverberation measurements and computations. The source was located at a 
depth of 48m with aU 16 elements receiving the reverberation. Seven different 
environmental models were used in the computation and were summarized in Table 1. 
The main differences among the seven models were the water column sound speed 
profile, interface roughness and volume perturbation. Of the seven profiles, one was 
chosen as the reference model with which the others were compared. Several analyses 
such as reverberation pressure levels, vertical correlation, peak correlation and spectral 
characteristics were performed using broadband signals. 

The transmission loss analysis provided a good preliminary prediction of what 
would be expected in the reverberation level analysis. It was found the three different 
sound speed profiles did not show much difference in the trans mi ssion loss. As for the 
variations in the interface roughness, the rougher interface produced less transmission 
loss, i.e. the signal was able to propagate further. This may be because the higher 
roughness values have caused the sediment layer between the water/bottom and the 
bottom/sub-bottom interfaces to be thinner and allow the acoustic energy to fateract more 
readily with the denser and faster volume below the bottom/sub-bottom interface. Thus, 
more forward propagation of the energy occurs. However, the higher volume 
perturbation increased the transmission loss significantly. This was probably due to the 
higher chances of producing a relatively slower bottom, thus, allowing the signal to 
penetrate more readily and not be reflected back into the water column. 

The variations in sound speed profiles of the water column showed little 
difference in the reverberation levels, correlation and spectral analyses. The variations in 
interface roughness showed that a rougher interface produced higher reverberation levels 
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and higher normahzed power in the power spectral analysis. This was consistent with he 
transmission loss analysis. In general, vertical correlation analysis relative to a source 
depth at 48m showed that the volume decorrelates much faster followed by the 
bottom/sub-bottom interface, then the water/bottom interface. It was also noted that the 
higher interface roughness caused the peak vertical coherence of the bottom/sub-bottom 
interface to be lost in the volume. The variations in volume perturbations showed that the 
higher sound speed perturbation produced lower reverberation levels and lower 
normahzed power in the power spectral analysis, which was consistent with the 
transmission loss result. There was no significant difference in the correlation analysis 
for the different volume perturbations, except when a large rms value was used. In the 
PRSD analysis, higher normalized power was generated by higher volume perturbation. 
There was no clear connection of this result with the transmission loss result at this time. 

With the conclusion of this thesis, recommendations for future work are : 

• To perform short-range statistical analysis to narrow down the relationship 
between signal stmcture and environmental stmcture. This is an attempt to 
interpret influence of multipath effects on long-range stmctures. 

• To perform vertical correlation analysis using a range of pulse lengths. This 
may help to distinguish interface from volume reverberation stmctures. 

• To incorporate environmental measurements from actual ASIAEX data and 
conduct prediction analysis. 

• To conduct data processing of ASIAEX data and perform data/model 
comparisons. 

• To incorporate rough sea surface scatter into the propagation model and 
investigate influence on various reverberation level predictions. 
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